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Abstract 

Ever since real- world two terminal memristor devices were created, de- 
tractors have derided them as not 'true' Chua memristors because there is 
no magnetic flux in the system. This report shows that the flux is present 
and equal to the magnetic flux of the oxygen vacancies. This shows that 
the HP memristor is a 'true' Chua memristor. The memristor description 
has been separated out into the memory function, which is the memrsi- 
tance described by Chua in 1971 and the conservation function, which is 
necessary to describe the whole HP memristor. As the memory function 
deals with the vacancy mobility and the measurable current is mostly elec- 
tronic, it is apparent that a description of the HP memristor must include 
two charge carriers. Prom analysis of Stan William's model in terms of 
memory and conservation functions, a direct relation to Chua's original 
memristance equation and an approximation of magnetic flux have both 
been identified as being present in the original model. With this, the 
phenomenological model used by experimentalists and the mathematical 
model beloved by theoreticians have been combined into one. 

1 Introduction 

There are four fundamental circuit properties which describe a circuit's opera- 
tion: the electrical potential difference, V, the electronic current, /, the magnetic 



flux, (yj and the charge, q. Three pairs of relationships were known in 1971 which 
defined the operation of the three (then) fundamental circuit elements: the re- 
sistor [R — V/I)\ the capacitor (C = q/y)] the inductor (/ = ^p/V). Analysed 
like this, there is a missing relationship that should relate charge and magnetic 
flux, and this was the relationship that Chua proposed [\\ would be satisfied by 
the discovery or creation of a memristor. From electrodynamics it is known that 
a moving charge has an associated magnetic field, so this idea seemed sound. 
Nonetheless, between 1971 and 2008, no one who had read Chua's work was able 
to make a memristor, perhaps because they were investigating magnetic mate- 
rials used for inductors, and the idea of a memristor languished in the drawer 
of theoretical curiosities. 

Over this time period, many experimentalists had reported 'anomalous' I- 
V curves with variations on a pinched hysteresis loop. As current is the time 
differential of charge (/ = dq/dt) and voltage is the time differential of magnetic 
fiux [V = dtp/dt), any linear relationship between charge and magnetic flux 
would manifest itself as a non-linear relationship between current and voltage, 
such as that which might describe a pinched hysteresis loop. It wasn't until 
2008 that Stan William's group at HP published a paper describing a working 
memristor [2 , complete with a phcnomcnological model for its operation and 
references to Chua's theoretical work. There were two problems however, the 
magnetic flux was missing from the model and the device did not have an 
appreciable magnetic field, as was expected. 

In that 2008 paper, Strukov et al stated that the missing magnetic flux did 
not matter, because all the definition of a memristor requires is the existence of a 
non-linear relationship between charge and flux to fit the mathematics. This pre- 
emptive explanation did not prevent some people [3] [1] from deriding the both 
the HP memristor as not a proper Chua memristor and the phcnomcnological 
model as a 'toy' model seemingly unrelated to Chua's equations. 

In this report I will demonstrate that the HP memristor does include mag- 
netic flux and using electrodynamic theory I will derive how the memristance 
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relates this flux to the charge flow. From there, I wih extend the description 
to cover the electronic current and voltage of the memristor and demonstrate 
that Stan William's 'toy' model actually includes an (admittedly accidental) 
approximation for the magnetic flux. 

2 Memristor behaviour. 

Chua asserted that a memristor would be a device governed by the following 
relationship [1]: 

dip = M{q)dq (1) 

where M is the memristance and is a function of q. If q doesn't change over 
time then the memristance is constant and Ohms law is sufficient to describe 
the device (ie it is a resistor). This is the theoretical model of a memristor. 

The HP memristor is a layer of titanium dioxide sandwiched between two 
electrodes and it works by being able to reversibly inter-convert between the 
high resistance, i?off, stoichiometric form Ti02 and the low resistance, Rom 
doped form, Ti02-i^, where x oxygen atoms (per mole) have been lost from the 
structure leaving a positive (p-type semiconductor) vacancies. Based on these 
quantities, the Strukov et al put forward the phenomenological model [2] of their 
memristor's memristance in a form similar to: 

M{q)^Ros{l-RonPq{t)) , 

where the actual titanium dioxide properties (ion mobility and physical di- 
mensions of the Ti02 layer) have been gathered into /?, the material parameter. 
NB. Strukov et al never expressed their model in this exact form. 

Memristive systems (Chua's generalisation of his ideas in 1976 |S]) have a 
hysteresis because at a given voltage there is more than one possible current, ie. 
there is more than one possible resistance. The value of the current measured 
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at a point must be related to the history of the device, if the memristive system 
is to have a memory. To make a memristor, rather than a memristive system, 
the device property that causes this change in resistance must be controUed by 
vohage - this is necessary for the memristor to be a two-terminal device (part 
of Chua's definition for a fundamental circuit element memristor). 

I'm going to approach the question of how the HP memristor can be a 
true Chua memristor backwards, by asking how a Chua memristor would work 
using knowledge of its behaviour gleaned from the Chua's equations. We're 
going to focus on the property in the memristor responsible for its memory, 
which is also a property that changes in response to voltage which I shall call 
the 'memory property'. I postulate that for there to be a memory in either a 
memristive system or a memristor, the memory property must be both separate 
from and slower to respond to a voltage change than the conducting electrons. 
This slower response time leads to the lag in current which gives rise to the 
hysteresis loop and explains the frequency dependence of memristance: if the 
voltage changes too fast for the memory property, it can't respond quick enough 
for a measurable change and the size of the hysteresis loop shrinks to a straight 
line (this is the Ohmic regime). The memory property has to respond to the 
voltage, which suggests that it either needs to be affected by the potential 
difference and therefore be charged or to undergo a structural change due to the 
electrical energy supplied. And, for the memristor to be of any real use, this 
change in memory property has to be (at least qualitatively) reversible, so the 
device can switch back and forth. 

There are several different possibilities for memory properties, such as charged 
ions in the PEO-PANI memristive system ^6 , or the concentration of spin elec- 
trons 7J in spintronic systems or the 'thermal' phase change that can be trig- 
gered by voltage in a VO2 thin film [5]. 

In the HP memristor the memory property is the charged oxygen vacancies 
in the Ti02-2:- This may be a strange way of thinking about hole dopants, as 
the hole is caused by the lack of an electron, however, there are three important 
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points about these hole dopants: 1, the conducting electrons move much faster 
than the dopants drift; 2, changing the dopant concentration in a volume of 
the memristor will change the resistance, conducting electrons won't do this; 
3, charged dopants will respond directly to voltage. Therefore, they obey the 
conditions presented above. 

If there is a second charge carrier in the memristor, then we have to ask 
the question of which charge is related to the magnetic flux (a question the 
HP group never thought to ask - they assume that q is the electronic charge). 
I am going to demonstrate that the relevant charge is the charge of the non- 
electronic charge carriers and therefore it is the relationship between this charge 
and its associated magnetic flux which defines the memristance. To be explicit, 
qy is the charge associated with the non-electronic charge carriers, this is the 
product of the formal charge on that type of charge carrier and the number of 
them present. Thus, although it is the eff'ect on the electronic current which 
is mcasiircd (and will be of use in real world devices), the electronic current is 
irrelevant in the actual process of memristance. 

This paper will calculate the magnetic flux of the non-electronic charge car- 
riers in a memristor and show that this fits Chua's definitions for a memristor. 
We will then demonstrate that this flux is present in both the real world HP 
memristor and in the model equations, by doing this, Chua's mathematical 
theory will be united with Strukov et al's phenomenological model. This new 
description will offer a new way of understanding the workings of real world 
devices, enabling better predictions of memristive materials. 

3 Method 

First, the method for calculating the magnetic flux for any memristor where 
non-electronic charge carriers are responsible for the memory property will be 
described, then the magnetic flux in the HP memristor will be derived. 

We start with the Biot-Savert law for the magnetic field associated with a 



5 



volume current. This is the most appropriate formulation of the Biot-Savert law 
because we are going to consider the magnetic flux just above the memristor 
surface where the memristor is best viewed as a 3-dimensional object. The 
Biot-Savert law comes from magnetostatics, a branch of theory that describes 
the magnetic effects due to constant currents, although our current will change, 
magnetostatics is still a valid approach because the changing current is far slower 
than that to which such theory is successfully applied (namely mains A.C. (50- 
60Hz)). 

The volume current, J, is given by 

~ Vol 

where qy is the charge in that volume due to the non-electronic charge car- 
riers, /i„ is the ion mobility of the non-electronic charge carriers and L is the 
average electric fleld causing the movement of charge. The magnetic fleld or 
magnetic flux density, B, at a point, p, associated with with this charge is given 
by: 

= ^ y -^dr (2) 

where /io is the permittivity of a vacuum, dJ and df are the unit vectors for 
J and r where r is the vector of length r from the volume infinitesimal dr to 
point p, given by r = {rj, Vyj, rj}. 

To get the magnetic flux through a surface associated with this fleld, if, we 
need to take the surface integral 



>fi = jB-dA 

where dA is the normal vector from the surface infinitesimal dA. This 
method should work for any memristor where the memory property is related 
to moving non-electronic charge carriers. 
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4 Example derivation for HP memristor 



In the HP memristor, the non-electronic charge carriers are the vacancies created 



volume of the device that contains flowing vacancies. This varies with time, but 
at an instant in time, t, the magnetic field is given by 



where w{t) is the position of the boundary between doped and undoped 
titanium dioxide at time t and < w < D. 

The X-axis is taken as the direction of current flow between the electrodes, 
the limit of which is D, the device thickness which is lOnm 2 . The y and z axes 
are in the plane of the electrodes, with the limits E and F, the length and width 
of the crossed electrode area and are both 50nm in the crossbar memristor. We 
take the terms of our integral in the coordinate system for inside the memristor, 
ie: r^, ry and r^. 

The integral is solved using the technique of integration by parts, taking 
the cross product in the numerator, Jxr, as dg^r^, r^, r^) and the denominator. 



in the Ti02 layer. The magnetic field integral in equation |3] is taken over the 




dr.j;drydrz . 
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Figure 1: The HP memristor. The shaded area is the doped low resistance tita- 
nium dioxide and the unshaded region is the stoichiometric high resistance tita- 
nium dioxide. The arrow indicates the direction of vacancy movement through 
the material, these move in the x direction. The boundary between the two 
materials is w. The limits of the titanium dioxide layer in the y and z directions 
is E and F. As the vacancies move to the right along the x axis, the magnetic 
B field associated with them curls around in an anti-clockwise direction (not 
shown). The surfaces which cross magnetic field lines are those in the x-y and 
x-z planes. It is these surfaces that the magnetic flux is calculated over. 
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Py and Pz contain only the dimensions of the memristor, so even if they are 
not analyticahy simple, they are easy to calculate numerically. As expected of 
a magnetic field, the divergence of the feild is zero. The curl is non-zero as the 
field curls around the current in an anti-clockwise direction. 

For a memristor which is close to being full with the maximum number of 
vacancies (ie the limit) the field at point p is given by 

B(p) {0, ~6.37qVxz,6.37qVxy}, 

where V is the applied voltage, p — {x,y, z} and x, y and z refer to a second 
set of coordinates which are located outside the memristor whose unit vectors 
are i, j and k 

As B is also known as the magnetic flux density, to calculate the magnetic 
flux we need to pick a surface to evaluate over. It makes sense to choose a 
surface that correlates to one of the surfaces of the device. Picking the surface 
just above the device {0<x<D,0<y<E, z = F^, we use the surface 
normal area infinitesimal, dA, which is given by dA = {0, 0, ij}. As is standard 
in electromagnetism, we integrate over the entire area. In this case, the limits 
of the surface are taken to be the dimensions of the device. 

The final value for ip over this surface is 

(fi = f^Lq^yXyPz . (4) 

where xyPz is the z component of the magnetic B field and this is generalised 
to any surface in the appendix |A.2[ We've now shown that there is a magnetic 
flux associated with the HP memristor's operation. 

^The volume current is constrained within the memristor and can be written in terms of 
coordinates inside the memristor. The magnetic feild (as caused by the volume current) can 
only exist outside the memristor and therefore can be written in terms of coordinates from 
outside the memristor. The two sets are labelled differently here to avoid confusion. If the 
distinction is not made between the two sets, then it's possible that the inside coordinates 
might be integrated over twice, which would be wrong. Perhaps confusingly, the limits are the 
similar. The inside coordinates have the limits; < < D; < ry < E; < < F\. The 
outside coordinates can go from — oo to oo but must avoid the volume within the memristor 

■^Actually z = F + dF so the surface is just above the memristor, avoiding any surface 
effects. 
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Putting in real-world values for the device characteristics for the HP mem- 
ristor gives a flux of 2.44x lO^^^Wb. In contrast, the magnetic flux associated 
with the conducting electrons through the same surface]^ is -4.07x lO^^^Wb. 
This is in the opposite direction and approximately 100 000 times bigger than 
the vacancies' magnetic flux. This may explain why the magnetic flux associated 
with memristor function has not been experimentally measured. 

5 Relation to other memristor models 
5.1 The memristance 

Chua's original formulation of memristance, equation[l] is also the instantaneous 
or chord resistance [SJ , which is defined as the gradient of the I-V curve at that 
point in time and this is the ratio between dip and dq. From this definition and 
equation [9] we can immediately write the memristance, M as: 

M{q(t)) = f^DELfi^PMt)) (5) 

where M is a function of q{t) because the movement of the boundary over 
time is proportional to the amount of vacancies created: Pz{w{t)) cx Pz{q{t)). 

This is the memristance as calculated from the top surface. We get the same 
numerical value of memristance if we calculate it from the other four surfaces 
which cut magnetic feild lines, but if we calculate it from a surface parallel 
to the magnetic field lines we get the answer of zero. This means that the 
memristance is not calculable from any surface perpendicular to the vacancy 
current, ie surfaces in the y-z plane. 

This is intriguing as it raises the possibility that memristance could be a 
tensor quantity, or at the very least, some care must be taken in calculating it. 
For the rest of this discussion, we will concentrate on the memristance through 

^To get the number of electrons, we've assumed that the Ti02 is acting hke a metal and 
every titanium atom is giving up a conducting electron. To get the number of oxygens that 
can be lost, we're assuming that maximum of 3% of available oxygen atoms 
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the surfaces which do cross magnetic field Unes (note that regardless of which 
of these surfaces we chose, we will get the same answer because the y and z 
components of the B field are equal). 

Equation [5] can be considered as three separate parts: 

1. universal constants: this term includes the effects of the permittiv- 
ity of a vacuum on memristance. It's inclusion in the equation clearly 
demonstrates that magnetism is involved in memristance. 

2. experimental constants: DEL, where DE is the surface the flux was cal- 
culated over, in this case the top of the device, and the electric field, and 
it doesn't change over the experiment. L is the root mean squared electric 
field strength if we are calculating memristance in an A.C. circuit. 

3. /3, the material parameter: iJ-yPz, this includes the physical dimensions of 
the device, but it will change throughout the experiment as a result of 
the moving boundary caused by the movement of vacancies in the device. 
This is the only term that contains variables. 

5.2 Memory Function + Conservation Function view of 
Memristance 

When measuring a memristor it is conventional to measure the electronic cur- 
rent, not the ionic current. As the electronic current is many times larger and 
faster than the movement of vacancies, we can ignore the vacancy contribution 
to the total flow of charge. What is therefore needed is the memristance as 
experienced by the conducting electrons. The component of that memristance 
which is directly due to the changing resistivity of the doped material is given 

by 

M,-{e-) = CM,+ {v+), 
where C is an experimentally determined parameter for the material. Since 
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the ion mobilities for the electrons, fif,- , and the vacancies are experimen- 
tally measured quantities, it is predicted that C = . 

In the HP system, M„+ , is the Chua memristance, this is the direct relation- 
ship between magnetic flux and charge and it only holds for that portion of the 
magnetic flux that is related to the charge carriers that hold the memory, ie the 
memory property of the device. The memory function is the Chua memristance 
in a form relevant to the conducting electrons and thus M^-(e^) the memory 
function of the HP memristor. 

Now we have a description of the memristance of the non-stoichiometric 
Ti02-a;, in terms of the memory property, we need to account for the rest of 
the device. 

In equation [5] w was allowed to change as a function of time, and this is 
exactly what happens in the HP memristor if the system is modelled as a moving 
boundary between two types of interconvertable materials. The model is not 
complete unless we include a description of the stoichiometric Ti02. This is 
necessary to ensure that space is conserved in the memristor model, therefore 
the function that does that will be called the conservation function, -Rcon- 

In the HP memristor, the conservation function is simply the resistance of 
the stoichiometric Ti02. From the geometry of the device and resistivity of 
stoichiometric Ti02, PoS, we can write the conservation function as 



_ {D ~ w)pos 

rtcon - 

The total resistance as experienced by the conducting electrons, i?tot is then 



R.o.it)^^-^^^+M.-i,e-) (6) 

Where Rtot{t) can be called a memristance because it is a resistance that 
changes with time due to the action of charge. This is not the Chua memristance, 
but does contain it. The Chua memristance is the memory function term where 
the resistivity of the material can change. The conservation function will also 



12 



change as a result of the charge (because w is dependent on q), this changes 
the effective size of undoped part of the memristor, but the resistivity of the 
material within that volume does not change. Essentially, the difference between 
the two functions is whether the intrinsic resistivity of the material can change, 
in the memory function it can, in the conservation function it can't. 

Plotting this equation with test values gives the pinched hysteresis loop in 
I-V space indicative of mcmristance, see figure |2] Separately, both conservation 
and memory functions are sufficient to give rise to a memristive I-V curve. A 
plot of the memory function by itself will give the memristance due to the change 
in size and resistivity of doped part. A plot of the conservation function by itself 
will show memristance due to the change in size of the undoped part. The two 
functions can have different magnitudes and fundamental frequencies. For an 
individual system the terms may be of the same order and interact or either one 
of them may be responsible for most of the measured memristive effect. That 
depends on the resistivity of the two materials, however the second case is more 
likely as appreciable memristance is easier to see if the two materials have very 
different resistivities. In the HP memristor, the ratio between the undoped and 
dope resistivity was 160 . 

We shall call i?tot {t) the total memristance and it is interesting as it is similar 
to Stan William's formulation of memristance. 

Rtot {t) expressed in eqn |6] shows that the HP memristor is a perfect Chua 
memristor because the entire equation can be written to be a function of a state 
variable, w [5.. I've now shown that the HP memristor contains magnetic flux, 
can be described using Chua's equations and is a true Chua memristor, I will 
now go on to evaluate the model presented in Strukov et al's paper. 

5.3 Relation to Strukov et al's model 

Strukov et al's model can be expressed in terms of conservation and memory 
functions: 
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Figure 2: Simulated I-V curve based on the memory-conservation model of the 
HP memristor 

Rtot{t) — [Rln + ^off — RonRotf][l3q] 

where the first bracket is the conservation function and the second is the 
memory function. (NB, the R^^ term was dropped in their paper [5] because it 
is many times smaller than the other terms in the description). The conservation 
function arises from the description of the memristor as two space-conserving 
variable resistors, see appendix |A.4[ exactly as in the model presented here. 

The memory function in this model is interesting. To understand why, we 
need to look at the role of the material parameter. Strukov et al never put 
their material properties into a single parameter in this way, I've done it all 
the way through this document for a very good reason. Despite Strukov et 
al's assertions, and many other scientists agreement, that HP model does not 
contain magnetic flux, it does. 

The material parameter, /3 has the units of the inverse of magnetic flux, ie. 
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Wb"'^. This makes sense when we examine what it is. 



£)2 



/3 is the ion mobihty divided by an area. We know that movement of charged 
particles has an associated magnetic fiux, and is Strukov et al's accidental 
approximation for ilj^ If the device was a wide as it was thick, the area, D^, 
would be one of the faces of the memristor that the magnetic flux field lines 
would cross, ie one of the faces memristance is calculable from. Compared to 
the correct form of the memristance, Strukov et al have missed out both the 
electric field and the magnetic permittivity. Essentially, if the magnetic flux 
were on the same order as the flux of charged particles, and we could simplify 
the geometry-related terms in the theory then /?~^ could be taken as an 
approximation to the magnetic flux. 

It is possible to find a term that is similar to Strukov et al's (3 term in the 
model presented in this paper. This is the material parameter, fiyPz, for the 
top surface of the memristor. is the dimensions parameter, a large term 
that involves only the three dimensions of the Ti02 layer (this is the geometry- 
related term referred to above). However, it is in units of and thus the 
D^^ part of the HP's /3 is an approximation for the dimensions parameter. It 
is possible that the Strukov et al's model works because the /3 term was close 
enough, D, E and F are all on the order of lO's of nm, therefore the error 
introduced between using instead of is easy to miss (for a full discussion 



of Strukov et al's derivation and their errors, see appendix A.4). Note, as /3 
does not vary in the HP model, the conservation function is responsible for all 
of the measured memristance. 

As an approximation for the magnetic flux, j3 still quantitatively wrong. For 
/3 to be a magnetic flux it must change as the device charges or discharges. This 
could be fixed by replacing D with w. Far more worrying is that (3^^ is 1x10^'* 



^Strukov et al never claimed that any part of their model contained the magnetic flux, 
^the magnetic flux is at 90°to the ion flow, so a cross-product should be involved and this 
is the cause of the complexity of the geometry-related terms 
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larger than the actual magnetic flux. If we look at the magnetic field that this 
approximation for the flux would require, we expect a magnetic feild vector of 



length ^ which is a field on the order of 1x10^'* Tesla (see appendix A.5). To 
put that in perspective, the largest pulsed magnetic feild that has been created 
was 2.8 kT and the average neutron star puts out a field of 1x10^ T. There's 
no way the HP memristor is more magnetic than 10 million neutron stars. 

If we take to be an approximation for the magnetic flux, then it's ap- 
parent that the memory function {/3q) in their model is ^, which is the inverse 
of the Chua memristance for that system and thus Strukov et al's equations are 
built on Chua's memristor equations. Therefore not only is the HP memristor a 
true Chua memristor, Strukov et al's model is, in form, a true Chua memristor 
model. 



6 Memristors as magnetic resistors 

There are four circuit properties that Chua originally used to predict the exis- 
tence of the memristor, ip, q, V and /. Because electromagnetism contains both 
electricity and magnetism (the difference between the two fields depends a lot 
on where you're standing and how fast you're moving) every circuit contains all 
four properties. 

In the memristor, we can use these four circuit properties to describe the 
system in a novel way. We can designate the charge associated with the mem- 
ristive magnetic flux as the magnetic charge ^magnetic- Note, I am not imply- 
ing that these charge carriers are in any way more magnetic than any other 
charge carrier (neither am I claiming that they are magnetic monopoles, al- 
though Umul [lOJ has suggested that electrons in spin-ice memristive material 
behave like magnetic monopoles) . Instead, I am conceptually separating out the 
charge responsible for the memory function of the memristor from the electrons. 
In the HP memristor, (Jmagnctic is Qv f\ 

®We can imagine systems where the relationship is not so simple, for example where the 
interaction of several ions is involved. 



16 



If we've designated the and to be 'magnetic' properties of the system, 
it makes sense (by symmetry) to designate the current / and the vohage V 
'electric' properties. In this way, if a resistor is a relationship between V and 
/, then the memristor can be viewed as a magnetic resistor. The main point 
of doing this is to underline that there are two separate systems at work: the 
'magnetic' system, which is the Chua memristance relation between charge and 
flux and the 'electric' system, which is the resulting relation between current 
and voltage. This does mean that to properly characterise a memristor, we must 
measure the flux response to the charge, because only by knowing that would 
you be able to predict the memristor's behaviour, measuring V-I is not enough 
(for a discussion on this point see [9]). 

This analysis does require that wc start thinking of magnetism in a novel 
way. Every moving charge, no matter how slow or seemingly insignificant has an 
associated magnetic flux. The discussion presented in this paper demonstrates 
that we can not simply ignore these charges or their associated magnetic fluxes, 
because the proper electromagnetic description of memristance requires them. 
To identify memristive materials, we need to think less like physicists and elec- 
tronic engineers who label magnetic materials only those with large responses 
to magnetic fields, and more like NMR spectroscopists who know that almost 
every material will respond in some way to a magnetic field. Thus, in this con- 
text, we should now expand the label of magnetic materials from the traditional 
materials to include semi-conductors and other memristive materials. 

7 Conclusion 

The magnetic flux has been found for the HP memristor, allowing Chua's the- 
oretical framework to be united with Strukov et al's phenomenological model. 
The standard model for a memristor has been expanded and separated out into 
the memory function and space-conservation parts of the theory, which will aid 
greater understanding of how memristive devices work. By underlining the im- 
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portance of a second slow ion current, it will now be easier to understand and 
design niemristor materials. It has been demonstrated that the HP meniristor 
is a true Chua memristor and magnetic flux was present as an approximation 
in the original model of this device. The analysis offered here for the first and 
perhaps archetypical memristor can be applied to other memristor systems. The 
concepts of memory function and conservation function can be carried across to 
other systems. 

A Appendix of Supplementary Information 

A.l Integrating by parts to get an expression for B. 

The actual integration by parts is: 



JJJ f {ra:,ry,r^) {-£--£r-£rg {ra:,ry,r,)] dr^dr^drz = f {r^,ry,r^) g {r^,ry,r^)- 




Which we set equal to the simpler shorthand: 




(7) 



To solve the integral in equation |3j which is 





we set 



1 



u = 



{rl+rl + r1)2 



and 



dv = Jxr = {0, 



}■ 



Vol ' Vol 



Thus 
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du 



and 



"""^ ' 2Vol ' 2Vol ^' 



Substituting for the right hand side of equation ([?]) and simphfying gives the 
following vector F 



where 



udv = F = {F^,Fy,F,} 



-pindef 

y ~ 



2Vol 



- {r.ry {rl {Py + rlf + ^ (2r2 + r^) + rl {2r^y + + 2r\ 



-r^ arctan 



' X' y 



rzJrl+rl + rl ^ 



and 



19 



_ HyLqxy j -pindef\w\E\F\ 

' ~ 2Vol 



o 

1^ 7^ 7^ 

pindef ^ U ^ 



+ry arctan 



{rl + r-f )^ + rt (rg + 2r^) + rg (2r^ + 5rgrf + 2rf ))) 



ri + ri 



We then substitute in the hmits for the integration, rx\Q , r^j^ and rz\o and 
replace Vol by wEF to simphfy Py where P^^ = x^-P^"'''=^|o lo lo'- Thus, 



1 (f2 (£;2 + F2)' + (2^2 ^ ^2) ^2 (2i;4 + 5S2^2 ^ 2F4) j j 



2M;i;i^ ((U;2 + ir2^ (£2 + J^2) (y,2 + £;2 + J^2)i j 

arctan 



and similarly 



E 



2(w;2+£;2_^^2)2 

1 (wF (e^ {E^+F^f + {E^ + 2F2) + {2E^ + bE^F^ + 2^*) ) ) 



2wEF 



^ f wF \ 

+E arctan , . 

\Ey/w'^ + E^ + F"^ J 



20 



It is not a very elegant solution. 

To get the solution for B, F is multiplied by the permittivity of a vacuum 
constants that was outside the integral in equation |3j ie ^ and 



B{p) = '^"'^^^'^ {0,~xzPy,xyP^}, 
an 



where the Vol occupied by the volume current has been gathered up into 
dimensional parameters Py and Pz so all the dimensions associated with the 
size of the doped part of the device are in one place. 



A. 2 Integration over the surface. 

The general form of (p is given by 



ip = —LqfiyijPk 
in 



(9) 



where i and j are the two chosen Cartesian coordinates of the surface you 
wish to integrate over and P^, refers to the geometric term for the magnetic field 
component for the other Cartesian axis. This is the part of the B field lines that 
is perpendicular to the chosen surface and therefore the only part that would 
contribute to the flux. The values are given in table [1] 



Device surface 


Area infinitesimal 


Integral 


Value for HP 
memristor 


Top 


dA = {0, 0, m 


f^top = Jo Ja B ■ dAdxdy 


3.186xl0-i5g 


Bottom 


dA = {0, 0, -zj} 


— "Ej — T) 

(^bottom = /„ /o B • dAdxdy 


-3.186xl0-i5g 


Front 


dA = {0,ik},0} 


'/'front = /(, B • dAdxdz 


-3.186xl0-i5g 


Back 


dA {0,-tk,0} 


<Pback = /„/(, B • dAdxdz 


a.isGxio-i^? 


Left 


dA = {jfc},0,0} 


</5icft = Jo /□ B • dAdydz 





Right 


dA = {-jfc},0,0} 


(fright = /n /n B • dAdydz 






Table 1: Table for the magnetic flux as calculated from the different possible 
surfaces of the memristor. The magnitude is the same for any of the four surfaces 
which cross the magnetic field, the sign changes dependent on the direction of 
the field. The end surfaces, labelled left and right, are perpendicular to the 
flow of vacancies, and this parallel to the magnetic field and therefore have a 
magnetic flux of 0. 
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Integrals are done over x, y or z which are coordinates outside the device. 
It's best to think of the surface as hovering an infinitesimal above the actual 
surface of the device, to avoid edge effects. Note, that the values of the non- 
zero (/s's will change as the device switches. This dependence is shown above by 
keeping g as a variable but the limits in the x are taken a D and are thus the 
limiting case, in actual fact they could be less. 

A. 3 Memristor theory Comparison 



A table to compare the differences between my theory and Strukov et al's model. 



Quantity 


HP Model 


My Model 


Memristance 


-Roff — PRoSRonQ 


M{q{t)) = ij^oom^p^ 


13 






M(^) 


Ml. 


1 


Memory function, [Mem] 




M{q) 


Conservation function 


[Rl, + Ros - RonRos][Mem] 




dd{t) 
dt 


II T - II ^ - II * 


T — ^ 



* Note, the equation iJ^v^ — fJ'vjj which is utilized by Strukov et al in 
their derivation is incorrect, because the current, i, is the electronic current 
and the movement of the boundary, '^^Jf^ , depends only on the drift velocity of 
vacancies, which is fiyL 



A. 4 Discussion of Strukov et al's phenomenological model 
derivation | [2] 

Strukov et al's 2008 Nature paper [2] announcing the manufacture of the Ti02/Ti02-, 
at HP included a set of equations which modelled this system as two variable 
resistors linked by a rule that ensured space conservation. The terseness of these 
equations, due to the necessity for brevity, obfuscate some assumptions made. 
For this reason, the derivation is discussed below. 

First, they define the position on the a;-axis of the boundary between the 
Ti02/Ti02-x as being equal to w, where w is restricted to varying between 
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and D, the length of the device, ie, < w < D. This makes sure that space 
is conserved, ie the total amount of Titanium dioxide can not change in an 
unphysical way. 

From their definition of w, they get the following equation for the memris- 
tance that relates I{t) and V{t): 

which they call equation 5 in their paper, and the fact that it is a distance 
divided by the total distance, D, means that the two resistances vary according 
to fractions of 1, and thus that space is conserved in the model. 

They then concern themselves with the change in w over time, They ex- 
pect that the boundary will move over time, dependent on the ion mobility(/i,,), 
which is the effective speed of an ion, in this case oxygen vacancies, in a unit 
electric field. The electric field is given by E in these equations, L in the deriva- 
tions in this paper. Thus, 

dw{t) _ 

which they do not state in the paper. 

Voltage, V , is equal to the electric fcild, measured between two points divided 
by the distance between those two points. If we look at the voltage across the 
two ends of the device (which is handily measurable in the lab), then we get, 
V = ED and so 

dw{t) _ V 
which is also not included in the paper. 



However, the next step is, and the following equation is given as equation 6, 

w{t 
dt 
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which introduces ^^^§7^ and follows on directly from equation 5 as a statement: 



dw{t) 



= IJ'V 



iR 

Id 



dt 



In the paper they state that they are assuming 'ohmic electrical conduc- 
tivity', namely that V = iR and so can be replaced with it. I have used the 
lower case i to indicate the current, mostly to show that it is not necessarily 
/. li V = iR, then this current must be the current that is related to the 
voltage which is measured between the two ends of the memristor, ie the elec- 
tronic current, Zg-. In the HP memristor there arc actually two currents at 
work, the electrical current, which is the electrons flowing through the device, 
and the current due to the oxygen vacancies, io+, whose movement changes the 
resistance. 

The second thing that is worrying about this equation is the resistance, R. 
What is it? There are two resistances in the system, Ron and -Ron- Why do the 
authors chose i?on? The answer is that they have decided that on the x axis 

corresponds to the side of the device where the doped, (i?on) material is. So 
the oxygen vacancy only passes through i?on material to get to the boundary. 
As most the potential drop is across the Ti02 material, and because i?off is 100 
times bigger than i?on, by approximating the voltage at the boundary as V they 
are introducing a small error of about 1 in 100 (as the voltage drop across Ti02 
is 100 times more than Ti02_x- 

To get the, now-famous, equation for memristance in terms of the charge 
transferred, q{t), they now integrate the above expression with respect to time: 




as j^i{t) = q{t), 



they write the answer as eqn 7, which is: 



w{t) = jl. 



i?, 



•on 



q{t). 



'V 



D 
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This is odd as i{t) hasn't been explicitly defined but from equation 6 where 
V is replaced with iR, this must be the electronic current. In that case, q{t) 
must be the amount of electronic charge that has passed through the device, 
Qg-. This is fine, except that Chua's formulation seems to suggest that the 
charge is transferred onto the device somehow. In fact it is, the Ti02 undergoes 
a chemical change which turns it to Ti02_a; as a result of something in the 
system and, over time, the amount of Ti02_a; increases. This is an increase of 
charge in the system as each vacancy is a positive charge. These charge carriers 
are what remembers the state of the memristor and that they remember it 
with zero power draw and thus can retain a state when the power is off. This 
argument strongly implies that the charge in the memristor equation should be 
that of the oxygen vacancies, not the electronic charge. 

The equation above is put into eqn 5 in the paper which gives: 



which because i?on is 100 times less than i?off, they disregard the first term 
and report the rest in their paper. 




and thus memristance, M, is equal to 




A. 5 B field associated with HP's V' 



If 



D2 



and 




B.dA. 
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and taken to be an approximation for the surface of the memristor xy, 
ie the area we integrate over then 



dA 



As dA = dxdy 



B 



dxdy dxdy dxdy /i„ 
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